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We present the proceedings of the lectures given at the 2008 Helmholtz International Sum- 
mer School Heavy Quark Physics at the Bogoliubov Laboratory of Theoretical Physics 
in Dubna. In two lectures we present recent results from the existing B-factories experi- 
ments, Belle and BaBar. The discussed topics include short phenomenological motivation, 
experimental methods and results on B meson oscillations, selected rare B meson decays 
(leptonic, b — » S7 and b — > s£ + £~), mixing and CP violation in the system of D° mesons, 
and leptonic decays of D s mesons. 

1 Introduction 

The lectures presented in this paper are a part of the B-factories lectures prepared in col- 
laboration with A.J. Bevan (also given in the proceedings of the school, p]). To obtain an 
approximate overview of recent results on flavour physics arising from Belle and BaBar both 
sets of presentations (each composed of two one-hour lectures) should be consulted. 

The lectures presented here include - beside the experimental methods and results - some 
short phenomenological sketches of motivation and/or interpretation of individual measure- 
ments. The author, being an experimentalist, should warn the reader that some examples of 
phenomenological interpretation are simplified and that serious theoretical treatment requires 
consultation of references given in the text. Examples are thus to be treated with a grain of 
salt; to quote the famous poet: "It is a curious fact that people are never so trivial as when 
they take themselves seriously." (O. Wilde, 1854 ~ 1900). 

A large majority of results presented in the lectures arise from the measurements performed 
with the two experiments taking data at the B-factories, e + e~ asymmetric colliders running 
at the centcr-of-mass (CM) energy y/s = mx(4s)C 2 0. T(4S*), a bb bound state with a mass 
just above the threshold for T(4S) — > BB decay, is a copious source of B meson pairs. Mesons 
are produced almost at rest in the CM system, but since the electron beam has an energy 
higher than the positron one, they are boosted and decay time dependent measurements of 
meson decays are thus possible. The Belle detector [2J, operating at the KEKB collider [3] 
in Tsukuba, Japan, has so far recorded an integrated luminosity of around 860 fb _1 , roughly 
corresponding to 950 x 10 6 pairs of B mesons 0. The BaBar detector g] at the PEP-II collider 
in Stanford, USA, has recorded around 550 fb _1 of data. 



1 Here and in the following we adopt a notation where rax represents a nominal mass value of particle X . If 
we refer to the reconstructed invariant mass of a system Y we use the notation m(Y). 
2 Both, B°B° and _B+ B~ pairs are produced, at approximately the same rate. 
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Beside the production of B meson pairs from the T(4S) other processes take place in e + e~ 
collisions at the given CM energy. For the subject of the lectures the most important is the 
continuum production of cc quark pairs, arising in e + e~ — > 7* — > cc. This is sketched in Fig. 
[U where the cross-section for hadron production in electron-positron collisions is plotted as a 
function of the CM collision energy. The cross-section for the production of cc pairs is larger 
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Figure 1: Cross-section for production of hadrons in e + e~ collisions as a function of y/s. The 
resonant production of T(45 f ) — > is represented as the light shaded area, and the continuum 
e + e~ — > 7* — > as the dark shaded area. 

than the one for the B meson production, at the integrated luminosity of KEKB it corresponds 
to around 1.1 x 10 9 produced pairs of charmed hadrons. 

In the course of the lectures we will mention also some related results from experiments other 
than B-factories, specifically the ones from the CDF-II experiment at Tevatron [S], recording 
data in pp collisions, and Cleo-c experiment [6] at the e + e~ collider CESR, running at the DD 
meson pair production threshold. All these experiments provide for a truly diverse experimental 
environment to study various aspects of heavy flavour physics. "We all live with the objective 
of being happy; our lives are all different and yet the same." (A. Frank, 1929 - 194-5). 



2 Lecture I 

"Never loose an opportunity of seeing anything beautiful, for beauty is God's handwriting." 
(R. W. Emerson, 1803-1882) 

2.1 B meson oscillations 

The mixing of neutral mesons, that is the transition of a neutral meson P° into its antiparticle 
and vice- versa, appears as a consequence of states of definite flavour (P°, P°) being a linear 
superposition of the eigenstates of an effective Hamiltonian (states of a simple exponential time 
evolution) Pi 2' 

\P h2 )=p\P°)± q \P ) . (1) 

For a thorough derivation of the equations describing the oscillations of mesons the reader is 
advised to follow [7]. 
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While the mass cigenstatcs have a simple time evolution, the time dependent decay rate of 
flavour cigenstatcs depends on the mixing parameters x and y, expressed in terms of the mass 
and width difference of Pi,2 as x ~ (mi — 7712)/? and y = (Ti — T2)/2T. T is the average decay 
width of the two mass eigenstates. The decay rate of a state initially produced as a P° is 

dT(P° f) 



dt 
<1 



\Atf + \±A f n cosh yt +(\A f \' 
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(2) 



In the above equation t is a dimcnsionlcss decay time, defined in terms of a proper decay time 
t' as t = t'T. The notation Af, At is used to represent instantaneous amplitudes for P° — > / 
and P° — > / decays. It is obvious from Eq. (|2|) that using the decay time distribution of 
experimentally accessible flavour eigenstates one can determine the mixing parameters x and y. 
Moreover, the effect of the mixing parameters on dT/dt depends on the chosen decay channel 
(Af, Af). The decay time distribution of an initially produced P° is obtained from Eq. @ by 
replacing Af <-> Af and q/p — > p/g. The decay time distributions for decays to conjugated final 
state / are obtained by a simple / — » / transformation. The above decay rates are illustrated 
in Fig. [2] for several values of x and y. 
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Figure 2: Probability for an initially produced meson P° to decay at time t as P° (full curve) 
or P° (dashed curve). Qualitatively the left example roughly corresponds to the case of B® 
mesons, the middle one to the case of B° mesons and the right one to the case of D° mesons. 
Note the logarithmic scale on the right plot. 

The neutral B meson pair^f] from T(45 l ) decays are produced in a quantum coherent state 
with the quantum numbers corresponding to that of the T(45). Before the coherence is dis- 
turbed by a decay of one of the mesons, the pair is always in a B° — B° state. The decay rates 
given above are valid only after the first of the two mesons decays. To be used in measurements 
of B° mesons produced from T(45), the decay time t in Eq. ([2]) should thus be changed to At, 
the difference between the decay times of the first and the second neutral B meson (and the 
exponential factor should include |Af| instead of t). 

The experimental method of measuring B° meson oscillation frequencjQ x relies on a similar 
method as the one used for measuring the CP violation pQ. However, instead of CP specific 

3 We will use notation B° for mesons, while for the strange B mesons we will use a strict B^ notation. 

4 Strictly speaking experiments in B system measure the mass difference between the two eigenstates, Am. 
However, since the dimcnsionlcss mixing parameter x = Am/F can be more directly compared for different 
meson species, we prefer to use this. Similarly as for the notation of B mesons, we use Am and x for the B® 
mesons and Am s and x a for B® mesons. In lecture II we will use xjj and to denote the corresponding mixing 
parameters in the D° system. 
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final states, flavour specific final states of B meson decays are used (like B° — > J/ipK , K*° — > 
A' + 7r~), which allow to determine the flavour of the decaying B meson. The method is sketched 
in Fig. [3] The measured At distribution deviates from Eq. |J2J) due to several reasons: usage of 
flavour specific final state (Af = Aj= 0), negligible decay width difference (y -C 1), probability 
of wrong flavour tagging (w) and finite accuracy in determination of At (resolution function 
i? S ig(At)). Taking into account these corrections, the final expected decay time distributions 
are 



dT(B° -> /) 

dAt 
dT(B° -> /) 

dAt 



e _ l A *l| J 4/| 2 [l + (1 - 2ti;) cos(arAt)] ® R sig (At) 
e" |At| | J 4 / | 2 [l-(l-2w)cos(xAt)] (g> R sig (At) , 



(3) 



where the ® sign denotes a convolution. The resolution function is composed as a convolution 
of several Gaussian functions [5]. The average accuracy of At determination is around 1.4 ps 
(the lifetime of B° mesons is 1.53 ps [9]). 



B or B 




fully reconstruct decay 
to flavour specific final state 
=> flavour at A t = 



K" 



determine time 
between decays 



tag flavour of B from charges 
of typical decay products 
=> flavour at A t 



Figure 3: Illustration of the method used to measure the B° oscillation frequency x. 



The most precise single measurement of x |10j uses several flavour specific final states to 
reconstruct the signal B° meson decays. Results are presented in Fig. @] (left) in form of the 
asymmetry 



dT(B° -> f)/dAt - dT(B a -> f)/dAt 
dT(B° -»■ f)/dAt + dT(B° -»■ f)/dAt 



(1 - 2w) cos x At (g> R sig (At) 



(4) 



The average value of existing Am measurements 1 1 1 j . expressed in terms of x = Am/T, is 
x = 0.776 ±0.008. 

Calculation of (B°\H e f f\B°) matrix element, visualized by the loop diagram of Fig. [4] 
(right), results in [T!2] 
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Figure 4: Left: Result of B° oscillation frequency measurement [TU] shown in the form of the 
asymmetry (fj|. Right: Loop diagram describing the B° — > B° transition (top) and the short 
distance contribution to D° — > D° (bottom) transition. 

The equation is written using a subscript q to emphasize that the same relation is also appro- 
priate for the system of B® mesons. Using the measured value of the oscillation frequency for 
B° mesons one can determine elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, if 
the QCD parameters i]b, Bb , /s are knowrjf). Due to their non-perturbative nature these 
quantities are difficult to estimate (usualy the lattice QCD calculations, LQCD, are exploited) 
and result in a large uncertianty of CKM elements determination. The constraints from the 
measured value of Am on parameters (p, ff) used to parametrize the CKM matrix [13j are 
shown in Fig. [5] (left) [14] . Since 2006 the oscillation frequency is measured also in the system 
of B® mesons, x s = 25.5 ± 0.6 [15]. In the ratio of Am s /Am the QCD uncertainties cancel 
to a large extent. The measured ratio Am s /Am is thus much more constraining than Am 
constraint alone (Fig. [5] (left)), and actually at the moment represents the most constraining 
measurement for the p among various flavour physics studies. 

2.2 Leptonic B meson decays 

Measurements of charged B meson leptonic decays are interesting for several reasons: theoreti- 
cally they are easier to interpret compared to semileptonic and hadronic decays, within the SM 
the measured rates can potentially yield the value of the least known CKM element V u b, and 
they are sensitive to possible contributions of processes beyond the SM. A Feynman diagram of 
an arbitrary pseudoscalar meson leptonic decay is shown in Fig. [5] (right). The QCD effects are 
described by a single parameter fp, the meson decay constant describing the overlap of the two 
quarks wave function. The leptonic decays of a pseudoscalar mesons are helicity suppressed, 
the expected ratios of decay widths are 1 : 4 x 10~ 3 : 10~ 7 for the t, fj, and e decays, respectively. 
Despite the problems due to at least two undetected neutrinos in the final state the decays to 
r leptons are the only decays observed so far. 

The method of measurement consist of fully (partially) reconstructing the accompanying 

5 Function Sq(x) is known. 
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Figure 5: Left: Constraints in the (p,fj) plane arising from various measurement [14j . The light 
shaded region denoted by Ato^ represents the constraint from the B° oscillation frequency, 
and the dark shaded region denoted by Ato^&Ato s the constraint from the ratio of B° and 
Bg oscillation frequencies. Right: Feynman diagram of a pseudoscalar meson leptonic decay. 
The QCD effects arc described by the decay constant fp. Beside the SM W + contribution also 
particles not included in the SM (like the charged Higgs boson) may contribute. 



B meson using a large number of hadronic (semileptonic) decay modes. After the particles 
assigned to the tagging meson are successfully identified one searches for one or three charged 
tracks originating from the r decay. Finally the energy in the electromagnetic calorimeter 
{Eecl) not assigned to the particles used in the previous reconstruction is examined. Signal 
decays with only neutrinos left in the final state are expected to peak at Eecl ~ 0. The 
i?ECL distribution of selected events in the measurement by Belle [TB] is shown in Fig. [5] (left). 
The excess of events above the expectation from MC simulation at low values of -Eecl is the 
signal for the B + — > t + v t decays. From the fit to the distribution the branching fraction 
Br(B + — ► t + v t ) = (1.65 igjf ±0 37) x 10~ 4 is obtained, where the main contribution to the 
systematic error arises from the uncertainties in the shape of the -Eecl signal and background 
distributions. A similar measurement performed by BaBar [17] yields Br(B + — > t + v t ) = 
(1.2 ± 0.4 ± 0.4) x 10 -4 , and iclusion of the Belle measurement using the hadronic tagging 
|18j results in the average of all measurements provided by Heavy Flavour Averaging Group. 
Br(B+ -> t+v t ) = (1.51 ± 0.33) x 10" 4 HQ 0. 
Calculation of T(B+ ->t+v t ) yields [2D] 

r(5+^rV) = ^|K6| 2 /imBm^(l-^|-) 2 [l-^|-tan 2 ^ 2 , (6) 



6 After the school an updated average of Belle and BaBar results appeared in |19| . Br(B+ — » t + u t ) 
(1.73 ±0.35) x 10" 4 . 
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Figure 6: Left: Distribution of the energy in the electromagnetic calorimeter for selected B + — > 
t + v t candidate events. The excess of data (points with error bars) above the expectation from 
simulation (histogram) is the signal. Right: Excluded region (shaded) in the (m#,tan/3) plane 
arising from the measurement of Br(B + — > t + u t ) [16] . 



where the last factor in brackets is a correction due to a possible contribution of the charged 
Higgs boson. The measured value is in agreement with the SM expectation (using LQCD 
prediction f B = (216 ± 22) MeV [2T], and \V ub \ = (3.9 ± 0.5) x 1(T3 [20]) and allows to put 
constraints on the parameters [mjj, tan (3) in the two Higgs doublet models (tojj is the charged 
Higgs boson mass and tan/3 is the ratio of the vacuum expectation values). The constraints 
arising from the Belle measurement are shown in Fig. [6] (right). 

2.3 b — > S7 decays 

Decays involving the b — > sj transition cannot occur at the tree level in SM. Such a flavor 
changing neutral current (FCNC) is only possible as a higher order process and is thus sensitive 
to possible contributions of New Physics (NP). Some possible diagrams, within and beyond 
the SM, are shown in Fig. [7] (left). At the parton level the photon energy in the CM frame is 
approximately half of the b quark mass. Also at the hadron level E 7 is sensitive to m,b, which 
is important for determination of \V u b\ and \V c b\ from semileptonic B decays. 

There are both, theoretical and experimental difficulties in the measurements. The former 
arise since in all experimental methods there is a lower cut-off applied to E y . To determine the 
branching fraction, for example, one has to extrapolate the partial rate for E 1 > E cut to the full 
energy region using models, which introduces theoretical uncertainties. On the experimental 
side the efforts arc being made to lower the cut-off, but this makes problems due to the huge 
backgrounds even more severe (see Fig. [7] (right)). The name of the game is thus to suppress 
the backgrounds to an acceptable level; "Your background and environment is with you for life. 
No question about that.'' (S. Connery, 1930). Methods of reconstruction may be divided into 
inclusive, semi-inclusive and exclusive ones. In an inclusive measurement only the photon is 
reconstructed. From the total E 1 distribution of events recorded at the T(4S) peak an analogous 
distribution of events, recorded 60 MeV below the peak is subtracted. The latter represents 
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Figure 7: Left: SM (top) and some NP (middle and bottom) contributions to the b — > sj 
process. Right: Simulated photon energy distributions from various processes. The smallest 
shaded region is the contribution of b — > sj (note the logarithmic scale). 



only the photons arising from e + e~ — > qq (continuum) events, and if the distribution is scaled 
according to the integrated luminosity of both samples, the remainder after the subtraction 
represents the energy distribution of photons from B meson decays. A search is made for 
photon pairs consistent with 7r° or r\ decays and such 7's are removed from the selected sample. 
The remaining background is estimated using simulated samples but normalized using data 
control samples. 

Result of such an inclusive method is shown in Fig. [S] (left) [55]. E 1 distribution peaks at 




Figure 8: Left: Raw CM system photon energy distribution for inclusively reconstructed b — > S7 
decays [22] . Right: Differential branching fraction of B — > X s j as a function of E 1 obtained in 
semi-inclusive measurement [23 . 

around half of the b quark mass and is consistent with zero above the kinematic limit for B — > 
K'y decays, confirming the correctness of the subtraction procedure. To determine the branching 
fraction and the correct shape of the energy distribution one has to apply a deconvolution 
method to the raw spectrum, correct it for the efficiency of reconstruction, subtract a simulated 
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contribution of b — > e?7 decays (~4%) and make the transformation to the B meson rest frame. 
The partial branching fraction in the interval E 7 > 1.7 GcV is found to be Br{b — > s-f) = 
(3.31 ± 0.19 ± 0.37 ± 0.01) x 10~ 4 . The last uncertainty is due to the boost from the CM to 
the B meson frame, and the largest systematic uncertainty arises from the normalization of 
backgrounds other than 7r° and r\. 

As an example of a semi-inclusive measurement we present the analysis of B — > X s j by 
BaBar collaboration [23], where X s represent a sum of various decay modes with K^, Ks, 
and 77 mesons in the final state. The photon energy is in this method calculated from the 
invariant mass of the hadronic system m(X s ) which results in a better resolution compared to 
the measured photon energy in the electromagnetic calorimeter. The background is suppressed 
using neural network for the rejection of continuum events and vetoes for 7's from ir° and 77 
mesons. To calculate the branching fraction for B — ► X s ~f the number of observed events must 
be corrected for the fraction of decays not taken into account in the reconstruction (25% at low 
m(X s ) due a to non-inclusion of Kl, and higher at higher masses). The resulting differential 
branching fraction for E 1 > 1.9 GeV is shown in Fig. [5] (right). The integral rate in the 
£ 7 > 1.9 GcV interval is found to be Br(b -> 37) = (3.27 ± 0.18 ±g;|g ±g;?^) x 10" 4 , where the 
last error is due to the QCD parameters affecting the efficiency. 

The measured branching fractions impose limits on possible contribution of charged Higgs 
boson. The world average of inclusive branching fraction is Br(b — > 37) = (3.52 ± 0.23 ± 
0.09) x 10" 4 HI]- The 95% C.L. limit following from [Ml is m H > 300 GcV for any value of 
tan p. In all measurements also the first and the second moment of the photon energy spectra 
are determined. These can be expressed with the same QCD parameters entering also the 
determination of \V u b\ and \V c b\ in inclusive semileptonic B decays. Details of a simultaneous 
fit performed to photon energy spectrum in b — > 57 and lepton momentum and hadronic mass 
spectra in semileptonic decays to determine the elements of CKM matrix is described in |llj . 

2.4 b -> s£ + £~ decays 

Decays involving the b — ► s£ + £~ parton process are another example of a FCNC transition. 
From that point of view they are interesting for the same reasons as the b — > S7 decays. Again, 
the inclusive decays are theoretically easier to interpret than the exclusive ones. Nevertheless, 
a lot of work has been done in identifying the observables in exclusive decays, especially B — > 
K*i + l~ , for which the theoretical uncertainties are small [35]. Feynman diagrams contributing 
to b — > s£ + £~ are shown in Fig. [9] (left). The differential decay rate dT/dq 2 , where q 2 is the 
invariant mass of the lepton pair, can be described in terms of effective Wilson coefficients 
Cf ff , C| ff and Cff , which include the perturbative part of the process and thus dependence on 
heavy masses of SM particles mw, mt, as well as on possible NP masses tonp- The absolute 
value of the coefficient C^ s can be constrained from the measured rate of b — ► S7 process, while 
in b — > s£ + £~ additional information (on sign of C^* 1 , as well as Cg ff , C^) can be obtained due 
to the interference of the two amplitudes shown in Fig. [9] (left). NP could change the values 
of the Wilson coefficients as well as add new operators causing the transition. In exclusive 
B — > K*£ + £~ decays the theoretical description includes beside the Wilson coefficients also the 
non-pcrturbative part expressed by the form factors which are predicted with an accuracy of 
around 30% |26j . This uncertainty is significantly reduced in some observables arising from 
the study of angular distributions, like the lepton forward-backward asymmetry {Afb) and the 
fraction of longitudinally polarized K*'s (Fl). 

BaBar performed a study of B — > K*£ + £~ decays in [57]. The reconstruction proceeds 
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Figure 9: Left: SM Fcynman diagrams for b — > s£ + £~. Right: Energy substituted B meson 
mass for reconstructed B — > K*£ + £~ decays. The distribution is shown for a part of the full 
sample with low q 2 [27\ . 



through identification of a lepton pair or e + e~) with invariant mass not in the range 

of charmonium states J/tp or ip(2S). The K* can be either charged or neutral, reconstructed 
through Ktt, Ktt° and Kgir final states. The signal can be seen in the energy substituted 
B meson mass, Mes — (Ecm/2) 2 — (J2iPi) 2 ' where Pi and Ecu are the B decay products 
momenta and e + e~ collision energy, respectively, calculated in the CM frame (Fig. [5] (right)). 
The background is composed of combinatorial one (described by reconstructing events with 
/i ± e =F lepton pairs), hadrons misidentified as muons (in fJ^/M^ channel) and peaking background 
from B — > Dtt decays where the charmed meson decays to K*tt (vetoed by requiring the 
invariant mass of K*ix not to be consistent with a D meson). 

For the reconstructed events the distribution of the kaon helicity angle in the rest frame of 
K* is investigated to obtain the fraction of longitudinally polarized K*'s (Fl). Fl value is then 
used as an input to the fit to the distribution of the angle between the lepton and K* in the t + £ ~ 
rest frame. This angle follows a 1 + 3Fl + (1 — F£) cos 2 9( + (8/3) Aps cos9g distribution. The 
Afb is the lepton forward-backward asymmetry which can be predicted in terms of the Wilson 
coefficients. In Fig. [TU] (left) measured Ap b is shown as a function of q 2 for the measurement 
by BaBar as well as the most recent measurement by Belle collaboration [35]. The measured 
values are compared to the SM prediction and the expectation for the Wilson coefficient Cf 
of reversed sign. In general the measurements seems to be shifted to larger asymmetry values 
than predicted. 

Similarly as for the b — > S7 decays, also for b — ► s£ + £~ semi-inclusive measurements have 
been performed by summing up various hadronic decay modes of the strange quark system (if 5 
or K ± with to 4 pions) accounting for around 70% of the total decay rate. The average of 
the branching fraction measurements is Br(B — * X s £ + £~) = (4.50±};of) x 10~ 6 [TT] . 

Various measurements of FCNC can be combined to put constraints on possible NP contri- 
bution to Wilson coefficients. Within a Minimal Flavour Violation scenario these constraints 
are presented in Fig. CE] (right) [29]. Measurements of Br(B -> X s ~f), Br(B -> X s £ + £~), 
Br(K — > ttvv) and Br(B s — > are used as the input. The combination of measurements 
is consistent with the SM {5Ci = 0) although there are large areas corresponding to non-SM 
contributions possible. 
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Figure 10: Left: Measured Afb as a function of q 2 in B — > K*i + H~ decays. The data points 
marked with triangles (dashed-dotted error bars) are from [27] and circular data points (full 
error bars) from [28] . The shaded regions represent q 2 intervals not included in the measurement 
{J/4>, il>(%S) regions). Full line represents the SM prediction and the dashed curve prediction 
with Cy S = — C 7 ' . Right: Constraints on NP contributions to Wilson coefficients arising 
from measurements of various FCNC processes [29]. The SM corresponds to 8Ci = 0. Light 
shaded areas represent the 68% C.L. and dark shaded the 95% C.L. region. 



3 Lecture II 

"Charm is. ..a way of getting the answer yes without having to ask any clear question." (A. 
Camus, 1913-1960) 

3.1 D meson oscillations 

The second lecture is devoted to results in physics of charmed hadrons from B-factorics. Charm 
physics in recent years gained in interest of both, experimental and theoretical physicists, mainly 
due to new interesting results from Belle and BaBar. Both experiments are not only factories of 
B mesons but also of charmed hadrons (see Section[T]and Fig. [I]). Contemporary charm physics 
has a twofold impact: as a ground of theory predictions tests, mainly tests of LQCD, and as 
a self-standing field of SM measurements and NP searches. An example of the first kind are 
the measurements of charmed meson decay constants, to be compared to LQCD calculations, 
to verify those and thus enable a more reliable estimates of the CKM matrix elements from 
the measurements in the B meson sector. The outstanding examples of the second group 
of measurements are recent observations of D° mixing and searches for the CP violation in 
processes involving charmed hadrons. 

Neutral D mesons are the only neutral meson system composed of up- like quarks. Hence 
a different contribution of virtual new particles than in the mixing of other neutral mesons is 
possible in the loops of diagrams describing the D° «-> D° transition (Fig. 2] (right)). However, 
the short distance contribution to the mixing rate, illustrated by the box diagram, is extremely 
small. The reason is the effective GIM suppression; calculation of the amplitude for this transi- 
tion reveals [30] that it is proportional to V* s V* d V u dV us (m 2 — m 2 )/m 2 . Hence the amplitude is 
doubly Cabibbo suppressed, and furthermore arises only as a consequence of SU(3) flavour sym- 
metry breaking. The resulting oscillation frequency defined in Sect. 12. H is \xd\ = 0(1O~ 5 ). This 
unobservable effect is hindered by long distance contribution to the transition amplitude, for ex- 
ample from states accessible to both, D° and D° (e.g. D° — > K + K~ — > D°). This contribution 
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is difficult to estimate. Current calculations [31] within the SM predict \xd\, \vd\ < 0(1O~ 2 ). 
The result illustrates the order of magnitude of the mixing parameters to be expected in the D° 
system (compare to measured values of x, x s given in Sect. 12. ip as well as the large theoretical 
uncertainty of the predictions. 

The time evolution of an initially produced D° meson follows Eq. ^ with the simplification 
due to \x D \, \vd\ < 1: 

dT( V /} = e-Mf + l^^A f f + 0(x%, vl) . (7) 

The time integrated rate for an initially produced D° meson to decay as a D°, Rm = (x 2 D + 
y 2 D )/2 ~ 10~ 4 , is small and represents the reason for a 31 years time span between the discovery 
of D° mesons and the experimental observation of the D° mixing. 

There are several methods and selection criteria common to various measurements of D° 
mixing. Tagging of the flavour of an initially produced D° meson is achieved by reconstruction 
of decays D* + — * D°irf or D*~ — > D°tt~. The charge of the characteristic low momentum pion 
7r s determines the tag. The energy released in the D* decay, q = m{D*) — m(D°) — , has 
a narrow peak for the signal events and thus helps in rejecting the combinatorial background. 
D° mesons produced in B decays have a different decay length distribution and kinematic 
properties than the mesons produced in fragmentation. In order to obtain a sample of neutral 
mesons with uniform properties one selects D* mesons with momentum above the kinematic 
limit for the B meson decays. The decay time is obtained from the reconstructed momentum 
and decay length of D° meson, and the latter is obtained from a common vertex of D decay 
products and an intersection point of D° momentum vector and the e + e~ interaction region. 

Methods of measuring the mixing parameters as well as sensitivities depend on specific 
final states chosen. The first to be described are decays to CP eigenstate fcp- In the limit 
of negligible CP symmetry violation (CPV, described in Section 13. 2p the mass eigenstates 
Di 2 coincide with the CP eigenstates (in case of no CPV q/p = 1, sec Eq. {!])). In decays 
D° — ► fcp only the mass eigenstate component of D° with the CP eigenvalue equal to the one 
of fcp contributes. By measuring the lifetime of D° in decays to fcp one thus determines the 
corresponding 1/Ti or I/T2. On the other hand, flavour specific final states like K~it + have a 
mixed CP symmetry. The measured value of the effective lifetime in these decays corresponds 
to a mixture of 1/Ti and l/IV The relation between the two lifetimes can be written as [35] 

t(£>°) 

r(fcp) = 7— , 8 

1 + VfVcp 

where r(fcp) and r(D°) are the lifetimes measured in D° — > fcp and D° — > K~ir + , respec- 
tively, rjf = ±1 denotes the CP eigenvalue of fcp- The relative difference of the lifetimes 
is described by the parameter yep- Expressed in terms of the mixing parameters, yep reads 
[32] yep = yo cose/) — (1/2)AmXd sin0, with Am and </> describing the CPV in mixing and in 
interference between mixing and decays, respectively. In case of no CPV, Am = (f> = and 
Vcp = UD- 

The measurement of yep by Belle [33] represents the first evidence of D° mixing Q Number 
of reconstructed decays to CP-even states K + K~ and 7r + 7r~ were 110 x 10 3 and 50 x 10 3 , with 



7 Publishcd simultaneously with the measurement of D° —* K^n decays by BaBar |34| which also gives 
evidence of the mixing. 
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purities of 98% and 92%, respectively. A simultaneous fit to the decay time distributions of 
KK, 7T7T and Ktt decays was performed with i/cp as a common free parameter. In order to 
perform a precision measurement of lifetime in each of the decay modes a special care should 
be devoted to a proper description of t resolution function in various data-taking periods (for 
details the reader is referred to the original publication). The t distributions and the result 
of the fit are presented in Fig. [TTJ The quality of the fit (x 2 /n.d.f. = 313/289) confirms an 
accurate description of the resolution effects. The measured value of yep is (1.31±0.32±0.25)% 




-2000 2000 4000 -2000 2000 4000 -2000 2000 4000 

t (fs) t (fs) t (fs) 



Figure 11: Decay time distributions of D° — > h + h~ (h = K, w) [33]. Hatched histogram is 
the contribution of background obtained from m(h + h~) sidebands. Full line is the result of a 
simultaneous fit to all three distributions with yep as a common free parameter. 



and the largest contribution to the systematic uncertainty arises due to a possible small detector 
induced bias in the decay time determination, yep deviates from the null value by more than 
three standard deviations including the systematic uncertainty. This evidence is confirmed by 
a similar measurement performed by the BaBar collaboration [35], finding yep = (1.24±0.30± 
0.13)%. 

Another possibility to look for the effect of mixing represent decays of initially produced 
D a, s to a wrong-sign final state K + ir~. While the more abundant D° decays lead to the 
K~tt + charge combination, the wrong-sign combination can be reached through doubly Cabibbo 
suppressed (DCS) decays or through a D° — ► D° mixing followed by a Cabibbo favored (CF) 
D° decay. In order to separate the mixing contribution from the DCS decays an analysis of the 
decay time distribution must be performed. The f-dependent decay rate, dT(D° — > K + tt~) oc 
[Rd + VRDy'p>t + (1/4) (x^ + y'D)^ 2 ] e ~ i consists of three terms corresponding to DCS term 
(Rd), mixing term {x% + y'£) and the interference between the two. Additional complication 
in the interpretation of the result arises since the decay rate depends on parameters x' and y' 
which are the mixing parameters rotated by a strong phase difference between the amplitudes 
of CF and DCS decays. In [33] BaBar collaboration fitted the t distribution of around 4000 
reconstructed wrong-sign decays. Result of the fit is presented in Fig. fTS] (left) in terms of the 
allowed region in (x l2 ,y r ) plane. While the central value is in the physically forbidden region 
(x 12 < 0) the no-mixing point ((x' 2 ,y') = (0, 0)) is excluded by a confidence level corresponding 
to 3.9 standard deviations. Numerically they find x' 2 = (-0.22 ± 0.33 ± 0.21) x 10~ 3 and 
y' = (0.97 ±0.44 ±0.31)%. 

The method which allows for a direct determination of both mixing parameters, xp, and 
ypn is the study of decays into self conjugated multi-body final states. Several intermediate 
resonances can contribute to such a final state. In the recent measurement by Belle |36) the 
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Figure 12: Left: Likelihood contours in (x' 2 ,y') plane arising from the measurement of D° — > 
K + tt~ decays [53] . The lines correspond to 1-5 a C.L., the dot represents the central value 
and the cross the no-mixing point. Right: 95% C.L. contour of (xd,Vd) as determined in the 
time-dependent Dalitz study of D° — > Kstt + tt~ [36]. Regions arising from the fit with neglected 
and allowed CPV, as well as neglecting or incorporating the systematic uncertainty are shown. 

Ksk + -k~ final state was analyzed, where contributions from CF decays (e.g. D° — > A'*~7r + ), 
DCS decays (e.g. D° — * K* + n~) and decays to CP eigenstates (e.g. D° — > p°Ks) are present. 
Individual contributions can be identified by analyzing the Dalitz distribution of the decay. 
Due to the interference among different types of decays it is possible to determine their relative 
phases (unlike in D° — > K + n~ decays where the relative phase between DCS and CF decays 
cannot be determined). And most importantly, since these types of intermediate states also 
exhibit a specific time evolution one can determine directly the mixing parameters xd and yn 
by studying the time evolution of the Dalitz distribution. The signal p.d.f. for a simultaneous 
fit to the Dalitz and decay-time distribution is 

|(A^+^-|^ (t))| 2 = |^(TO 2 _,m^)[e- 4Al * + e- 4A2t ] + ^4(m 2 „ , m 2 ,) [e^ 1 * - e^ 2 '] | 2 , (9) 

composed of an instantaneous amplitude for D° decay, A{m_,m~,), and an amplitude for 
the 5° decay, ^(m 2 .,™^), arising due to a possibility of mixing. They both depend on the 
Dalitz variables to 2 = m 2 (Ksir~) and m 2 ^ = m 2 (Ks^ + )- The dependence on the mixing 
parameters is hidden in A1.2 = TOi,2 — £Ti 2/2. If CPV is neglected the amplitude for D° 
tagged decays is \(K s ir + ir-\D (t))\ 2 (m 2 + , ml, t) = \{K S TT + ^\D Q {t))\ 2 {m 2 _,m 2 + ,t). As in the 
case of B° oscillation measurements the p.d.f. of Eq. (J9j) must be corrected to include the finite 
resolution on the decay time. 

Based on ~ 500 x 10 3 reconstructed decays with a purity of 95% Belle obtained a good 
description of the Dalitz distribution using 18 different resonant intermediate states and a non- 
resonant contribution. A simultaneous fit to m 2 _, m\ and t yielded mixing parameters Xr> = 
(0.80 ± 0.29±g;^)%, y D = (0.33 ± 0.24±g;^)%. This represents by far the most constraining 
determination of xjj up to date. Contour of allowed (xD,yo) values at 95% C.L. is shown in 
Fig. QJ] (right). 
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3.2 CP violation in the system of neutral D mesons 

A general, easy to reach expectation is that possible CPV in processes involving charmed 
hadrons must be small within the SM. This arises due to the fact that such processes involve the 
first two generations of quarks for which the elements of the CKM matrix are almost completely 
real. Typical CKM factor entering both the short distance box diagram as well as the decays to 
real states accessible to both, D° and D°, is V* S V US - Using CKM matrix unitarity this can be 
expressed as ~V* d V u d[l + (V* b V u b) / (V* d V u d)\- Considering the small absolute value of the second 
term one can see that axg(V* s V U s) ~ $*((V* b V u b) / {VcdYud)) ~ 7x 10 -4 . This is the typical value 
of the weak phase in charmed hadron processes which determines the size of the CPV effects 
in the SM. For example, CPV asymmetries like Ar discussed below, are typically of the order 
of XdSuk/), where is the weak phase considered, and hence A? ~ O(10~ 5 ). Deviation of 
\q/p\ value from unity, which also represents the CP violation, is roughly expected to be of the 
order of sin</> ~ 10 -3 . These values are all below the current experimental sensitivity and any 
positive experimental signature would be a clear sign of some contribution beyond the SM. 

All three distinct types of CP violation, CPV in decays, in mixing and in the interference 
between decays with and without mixing (see lectures by A.J. Bevan) can in principle be present 
in the D° system. They are parameterized by Ajj, Am and <j) according to 

4 1 = 1 + 4?; {Ad + 0, CPV in decay) 

1-1 = i + 4r : ( Am * °> CPV in mixin §) 

P 2 

<M/ _ (l + i M /2)^^_ J/} 



pA f l + A D /2 



e W-6f). (0 ^ 0, CPV in interference) . (10) 



In the above equations \JRd is the ratio of amplitude magnitudes \Af/Af\ and Sf is the strong 
phase difference between the two. 

In all mentioned mixing parameter measurements also a search for possible CPV has been 
performed 0. 

In decays to CP eigenstates (D° — * KK, tttt) one measures lifetimes separately for D° and 
D° tagged events. A measurable asymmetry 

A _ r(D° ^ fcp) - r(D° -> f CP ) ( 
r T(D°^fcp)+T(D°^f C p) [ ' 

is related to the mixing and CPV parameters as [35] Ay = (1/2) AmUd cos0 — xosmcj) and 
equals zero in the case of no CPV. The measured values by Belle [33] and Babar [35] are 
A T = (0.01 ± 0.30 ± 0.15)% and A T = (0.26 ± 0.36 ± 0.08)%, respectively. Hence there is no 
sign of the CP violation at the sensitivity level of around 0.3%. 

In D° — * K + ir~ decays the decay time distribution is also fitted separately for D° mesons 
and their anti-particles. There are six observables, x /=t , y /=t , R^, where the ± superscripts 
denote the observables for D° and D° subsamples. They are related to the parameters of Eq. 
m^yAn = (R+-R D )/(R++R- D ),R M = (x' ±2 +y' ±2 )/2 andA M = (R+ - R M ) / (R+ + R M ) . 



8 Note that both, mixing and CPV searches can also be performed using the time-integrated quantities, 
for example the rate of wrong-sign semileptonic decays D° — ► D° — > t~ K + v [37| or the CP asymmetry 

(r(D° -» /) - r(D° - /))/(r(r>o - /) + r(£>° /)) Eg. 
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From such fits the results of the search for CPV in mixing and in decay arc Am = 0.1 ± 
2.9, A D = (-2.1 ±5.4)% [34] or Am =0.67±1.2, A d = (-2.3 ±4.7)% [39] (errors here include 
statistical and systematic uncertainties). There is no hint of a direct CPV at the level of 5%. 

In the i-dependent analysis of D° — > Kstttt Dalitz distribution the possibility of CPV is 
included by additional two free parameters in the fit, Am and cf>. Also the direct CPV can 
be checked by allowing the contributions of various intermediate states to be different for D° 
and D Dalitz distributions. The latter was not observed within the statistical uncertainties. 
Parameters of CPV in mixing and interference arc found to be Am = —0.26 ± 0.60 ± 0.20 and 
= (—0.24 ± 0.31 ± 0.09) rad. The contours of (ijj, yu) arising from the fit allowing for the 
CPV are presented in Fig. [12] (right). 

3.3 Average of D° mixing parameters 

To make conclusions arising from a variety of results on D° mixing and CPV searches the 
Heavy Flavour Averaging Group performs an average of various measurements including cor- 
relations among the measured variables [llj . An illustration of (xo,Hd) constraints imposed 
by individual measurements is shown in Fig. [J3] (left). World average of mixing and CPV 
parameters for the D° system is presented in Tab. [T] The results arc presented graphically 




Figure 13: Left: Illustration of constraints on {xdiDd) values arising from various measure- 
ments. Middle: Probability contours for (xd,Ud) corresponding to 1-5 a C.L. from the average 
of measurements QT|. Right: Probability contours for CPV parameters (\q/p\,(j>) corresponding 
to 1-5 a C.L. from the average of measurements [TT] . 



Parameter 


Value 


Parameter 


Value 


x D 


(1.00±£i)% 


A D 


(-2.1 ±2.4)% 


VD 


(o.76±£i5)% 


Wp\ 


0.86±^ 


Rd 


(0.336 ±0.009)% 





o ooi 7.6" 
— o.o nij 2° 



Table 1: Average of D° mixing (left) and CPV (right) parameters [TT] . 

in Figs. 1131 (middle) and 1131 fright) as contours in (xD,y£,) and (\q/p\,<p) planes. The mixing 
phenomena in the neutral D meson system is firmly established, with the mixing parameters x d 
and no of the order of 1%. The oscillation frequency can be compared to the values for other 
neutral meson systems, x w 0.8 (B°), xk ~ 1 (K°) and x s w 25 (B®). Since both parameters, 
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Xd and t/d, appear to be positive, it seems that the CP-even state of neutral charmed mesons 
is shorter-lived (like in the K° system) and also heavier (unlike in the K° system). At the 
moment there is no sign of CPV in the D° system, at the level of one standard deviation of 
the world average results. 

Results in the D° mixing impose some stringent constraints on the parameters of various 
NP models [3D]. As an example we quote the i?-parity violating Supersymmetry models, where 
an enhancement of xjj could arise from an exchange of down-like squarks or sleptons in the 
loop of the box diagram. The exclusion region of possible values of the squark mass and 
i?-parity violating coupling constants for various upper limits on xd is presented in Fig. fH] 
(left). Planned Super B-factory, which would accumulate data corresponding to an integrated 




Figure 14: Left: Constraints on the values of squark mass and i?-parity violating coupling 
constants arising from xd < 1% (hatched region is excluded). Dashed lines represent boundaries 
of the exclusion region for xd < 1.5%, 0.8%, 0.5% and 0.3% [J0J. Right: Sketch of a method 
to measure Br(Df — > n + v) [41]. Full lines represent particles detected in the detector or 
exclusively reconstructed, and dashed lines particles reconstructed in the recoil (from the known 
momenta of incident beams and detected particles). 

luminosity of 50 ab _1 (compared to the current 0.8 ab _1 at KEKB), would of course yield 
results on D° mixing and CPV of much better precision. The extrapolated accuracies are 
ct(xd) - 0.1%, a(y D ) ~ 0.06%, a(\q/p\) ~ 0.05 and a{4>) ~ 3°. This would allow to severely 
constrain relations among parameters of various NP parameters and to search for possible CPV 
phenomena in the region where a large number of these models predict an observable effect. 
However, one should not forget the words "Prediction is very difficult, especially of the future. " 
(N. Bohr, 1885 - 1962). 

3.4 D s leptonic decays 

Charmed mesons leptonic decays are analogous to the leptonic decays of B mesons (Sect. I2.2() . 
By measuring the rate of such decays one would hope to determine the decay constant of the 
corresponding meson, sec Eq. ([6]), and by that test the predictions of LQCD. Both Belle and 
BaBar performed measurements of Br(D~^ — > ti + v). Clco-c collaboration measured decays 
Df — > tv as well. 

In [44] Belle measured the absolute branching fraction of — * n + v using a method illus- 
trated in Fig. [14] (right). Events of the type e + e~ — > D*D ± K ±fi X are used, where X can be 
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any number of additional pions from fragmentation, and up to one photon. An event is divided 
into a tag side, where a full reconstruction of D and a primary K meson is performed, and a 
signal side where the decay chain D* + — > Df-f, Df — > is searched for. Tag side charged 

and neutral D mesons are reconstructed in D — > Krm decays. For all possible combinations 
of particles in X, the signal side D* + meson is identified by reconstruction of the recoil mass 
m lec (DKX), using the known beam momentum and four-momentum conservation. The recoil 
mass m roc (y) is calculated as the magnitude of the four-momentum Pbeams — Vy ■ The next 
step in the event reconstruction is a search for a photon for which the recoil mass m rcc (DKXj) 
is consistent with the nominal mass of . The sample of Df mesons reconstructed using 
this procedure represent an inclusive sample of decays, among which the leptonic decays are 
searched for. If an identified muon is found among the tracks so far not used in the reconstruc- 
tion, the square of the recoil mass m% cc (DKX^n) is calculated. For signal decays this mass 
corresponds to the mass of the final state neutrino and hence peaks at zero. 

Final distribution of m^ cc (DKX^^,) is shown in Fig. [15] (left) where a clear signal of leptonic 
decays can be seen. Majority of background can be described using reconstructed — > e + v 
decays where due to the helicity suppression no signal is expected. Number of reconstructed 
signal decays is found to be N(Df — > n + v) = 169 ± 16 ± 8. Comparing to the number of 
inclusively reconstructed Df decays and correcting for the efficiency of muon reconstruction 
one obtains the branching fraction Br(Df — > ^ + v) = (6.44 ± 0.76 ± 0.56) x 10~ 3 . The largest 
contribution to the systematic uncertainty arises from a limited number of simulated decays 
used to describe the shape of the signal distribution. Using Eq. © (without the factor arising 
from the charged Higgs contribution) and the value of |V CS | as determined in a global fit to the 
CKM elements applying the unitarity of the matrix [9] , one determines the value of D s meson 
decay constant, f Ds = (275 ± 16 ± 12) MeV. 

BaBar [32] used a somewhat different approach by measuring the yield of -D+ — ► [i + v 
relative to the Df — > 07r decays. The branching fraction determined in this way is a relative 
measurement normalized to the Br(Df — ► <fiir). While this method enables a larger statistics 
of the reconstructed sample it suffers from a hard-to-estimate systematic uncertainty in the 
normalization mode (4>tt state is actually an intermediate state of the K + K~n + final state 
and can be influenced by the interference among various intermediate states). The neutrino 
momentum is determined from the missing momentum in an event. The resolution is improved 
by constraining the v and the reconstructed muon momentum to yield the nominal mass of D a 
meson. Fig. [15] (right) shows the distribution of the reconstructed mass difference between the 
D* and D s meson, where the signal of leptonic decays consist of 489 ± 55 events (the error is 
statistical only). Calculation of the decay constant yields a value of fn s = (283 ± 17 ± 7 ± 
14) MeV, where the last error is due to the uncertainty of Br(Df — ^ tfiir). 

How do the measured values compare to the LQCD calculations? The average of absolute 
measurements (beside the described Belle measurement these include measurements by Cleo-c 
collaboration in muon and tau decay modes [44]) is fu t ~ (274 ± 10) MeV [20]. The recent 
LQCD result exhibits a huge improvement in the accuracy compared to previous determinations 
of the D s meson decay constant: fo s = (241 ± 3) MeV [43]. The discrepancy between the two 
values is more than 3 standard deviations. While the fact that for the D + decay constant 
experimental results confirm the calculation (albeit within larger errors) may point to some 
intervention of NP [45] one should probably wait for a) confirmation of the LQCD estimate 
(and especially its uncertainty) and b) more accurate experimental measurements before making 



18 



HQP08 



any conclusion 0. 




0.05 



0.1 



0.15 



0.2 



0.25 




0.3 o.:- 

AM (GcV/c 3 ) 



0.35 



ML(DKXin)/GeVV 



Figure 15: Left: Distribution of m^ ec (DKXjfj,) for Df — > n + v candidate events [41]. The peak 
at zero corresponds to the signal decays. Right: m(D*) — m(D s ) distribution for Z?+ — > n + v 
candidate events [42]. Signal peaks at the nominal mass difference. 



4 Summary 

Although in the lectures we were able to present only a small fraction of exciting physics results 
that arose from the B- factories over almost a decade of operation, we hope the selected examples 
demonstrate the following: 

• B-factories have successfully performed precision measurements in identification of SM 
processes and determination of SM parameters, as well as a complement to direct NP 
searches that are soon to be started at the LHC; 

• experimental tests in general confirm predictions of the SM, although several hints of 
discrepancies at the level of 3 standard deviations exist; 

• _B-factorics have outreached their program as foreseen at the startup. 
Specifically related to the presented measurements one should note 

• B oscillations in conjunction with a breakthrough in B s oscillations confirm the SM to a 
high accuracy; 

• leptonic and radiative B meson decays constrain possible contribution of NP but large 
room for improvement remains for the Super S-factory; 

• important results in charm physics complement the results in the B meson sector; 

• measurements of D mixing and search for the CPV represent another achieved milestone 
in particle physics, more precise measurements and theoretical predictions are needed; 

• D s leptonic decays may test predictions of LQCD once the results are confirmed. 

9 For D+ decays the contribution is proportional to (m 2 D /m? H )(m £ j/m c ) tan 2 f} while for decays it is 
proportional to (m 2 -, /m 2 H )(m a /m c ) tan 2 /? 1461 . See also Eq. (0, where for £?+ decays, due to mj 3> m«, the 
correction is simply (?n 2 B /m 2 H ) tan 2 (3. 
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